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At very high densities, electrons react with protons to form neutron rich matter. This material is 
central to many fundamental questions in nuclear physics and astrophysics. Moreover, neutron 



£^ rich matter is being studied with an extraordinary variety of new tools such as the Facility for Rare 



Isotope Beams (FRIB) and the Laser Interferometer Gravitational Wave Observatory (LIGO). We 
describe the Lead Radius Experiment (PREX) that uses parity violating electron scattering to 

i 1 measure the neutron radius of 208 Pb. This has important implications for neutron stars and their 

crusts. We discuss X-ray observations of neutron star radii. These also have important impli- 
i— h cations for neutron rich matter. Gravitational waves (GW) open a new window on neutron rich 

^ matter. They come from sources such as neutron star mergers, rotating neutron star mountains, 

and collective r-mode oscillations. Using large scale molecular dynamics simulations, we find 
neutron star crust to be very strong. It can support mountains on rotating neutron stars large 
enough to generate detectable gravitational waves. Finally, neutrinos from core collapse super- 
novae (SN) provide another, qualitatively different probe of neutron rich matter. Neutrinos escape 
from the surface of last scattering known as the neutrino-sphere. This is a low density warm gas 
^sO of neutron rich matter. Neutrino-sphere conditions can be simulated in the laboratory with heavy 

CNJ ion collisions. Observations of neutrinos can probe nucleosyntheses in SN. We believe that comb- 

■ 

ing astronomical observations using photons, GW, and neutrinos, with laboratory experiments on 
nuclei, heavy ion collisions, and radioactive beams will fundamentally advance our knowledge 
of compact objects in the heavens, the dense phases of QCD, the origin of the elements, and of 
neutron rich matter. 
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1. Introduction 

Multi-messenger astronomy observes matter under extreme conditions. In this paper we de- 
scribe how electromagnetic, gravitational wave, and neutrino astronomy, along with laboratory 
experiments, provide complimentary information on neutron rich matter. Compress almost any- 
thing to very high densities and electrons react with protons to form neutron rich matter. This 
material is at the heart of many fundamental questions in Nuclear Physics and Astrophysics. 

• What are the high density phases of QCD? 

• Where did the chemical elements come from? 

• What is the structure of many compact and energetic objects in the heavens, and what deter- 
mines their electromagnetic, neutrino, and gravitational-wave radiations? 

Furthermore, neutron rich matter is being studied with an extraordinary variety of new tools such 
as the Facility for Rare Isotope Beams (FRIB), a heavy ion accelerator to be built at Michigan State 
University [1], and the Laser Interferometer Gravitational Wave Observatory (LIGO) [2]. Indeed 
there are many, qualitatively different, probes of neutron rich matter including precision laboratory 
measurements on stable nuclei and experiments with neutron rich radioactive beams. While as- 
trophysical observations probe neutron rich matter with electromagnetic radiation, neutrinos, and 
gravitational waves. In this paper we give brief examples of how neutron rich matter is being 
studied with these extraordinarily different probes. 

We are interested in neutron rich matter over a tremendous range of densities and temperatures 
were it can be a gas, a liquid, a solid, a plasma, a liquid crystal, a superconductor, a superfluid, 
a color superconductor, etc. Neutron rich matter is a remarkably versatile material. The liquid 
crystal phases are known as nuclear pasta and arise because of coulomb frustration [3, 4]. Pasta 
is expected at the base of the crust in a neutron star and can involve complex shapes such as long 
rods ("spaghetti") or flat plates ("lasagna"). Neutrinos in core collapse supernovae may scatter 
coherently from these shapes (neutrino pasta scattering) because the shapes have sizes comparable 
to the neutrino wavelength [5]. An example of nuclear pasta is shown in Fig. 1. 

At very high densities, neutron rich matter should form exotic quark and gluon phases. Asymp- 
totic freedom in QCD implies that at high densities, quarks will be nearly free. Under these condi- 
tions, attractive gluon exchange interactions should pair the quarks into a color superconductor [7]. 
Multi-messenger astronomy is, at present, the only way to observe cold dense matter and search 
for a color superconductor. It can not be created in the laboratory. Although heavy ion collisions 
can produce high densities, there is no way to get the entropy out and produce low temperatures. 
Thus hot dense matter can be studied in the laboratory, but not cold dense matter. Indeed a very 
interesting strongly interacting quark gluon plasma has been observed at the Relativistic Heavy Ion 
Collider (RHIC) [8]. This very hot material forms a nearly perfect fluid with a low shear viscosity. 

In this paper we focus on some of the simpler gas, solid, and liquid phases of neutron rich 
matter. In Section 2 we describe a precision laboratory experiment called PREX to measure the 
neutron radius of 208 Pb. Nuclei are liquid drops, so PREX and many other laboratory experiments 
probe the liquid phase of neutron rich matter. In astrophysics, electromagnetic, gravitational wave, 
and neutrino probes can observe different phases of neutron rich matter because the probes have 
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Figure 1: Surfaces of proton density for a pasta configuration of neutron rich matter at a baryon density of 
0.05 fm~ 3 . This is from a semiclassical molecular dynamics simulation with 100,000 nucleons [6]. 

very different mean free paths. In Section 3 we describe electromagnetic observations of neutron 
star radii. In Section 4 we discuss gravitational waves from neutron star mergers that are produced 
by the energetic motions of dense liquid phase neutron rich matter. In addition, continuous gravita- 
tional waves can be produced by "mountains" of solid neutron rich matter on rapidly rotating stars. 
In Section 5 we discuss neutrinos from core collapse supernovae. The neutrinos are emitted from 
a low density warm gas of neutron rich matter. We conclude in Section 6. 

2. Laboratory probes of neutron rich matter 

Neutron rich matter can be studied in the laboratory. Hot and or dense matter can be formed in 
heavy ion collisions, while more neutron rich conditions can be accessed with radioactive beams. 
In addition precise experiments are possible on stable neutron rich nuclei. We give one example, the 
Lead Radius Experiment (PREX) [9] accurately measures the neutron radius in 208 Pb with parity 
violating electron scattering [10]. This has many implications for nuclear structure, astrophysics, 
atomic parity violation, and low energy tests of the standard model. 

2.1 Introduction to neutron densities and neutron radii 

Nuclear charge densities have been accurately measured with electron scattering and have 
become our picture of the atomic nucleus, see for example ref. [11]. These measurements have 
had an enormous impact. In contrast, our knowledge of neutron densities comes primarily from 
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hadron scattering experiments involving for example pions [12], protons [13, 14, 15], or antiprotons 
[16, 17]. See also ref. [18] for a beautiful measurement of the dipole polarizability. However, the 
interpretation of hadron scattering experiments is model dependent because of uncertainties in the 
strong interactions. 

Parity violating electron scattering provides a model independent probe of neutron densities 
that is free from most strong interaction uncertainties. This is because the weak charge of a neutron 
is much larger than that of a proton [19]. Therefore the Z° boson, that carries the weak force, 
couples primarily to neutrons. In Born approximation, the parity violating asymmetry A pv , the 
fractional difference in cross sections for positive and negative helicity electrons, is proportional to 
the weak form factor. This is very close to the Fourier transform of the neutron density. Therefore 
the neutron density can be extracted from an electro-weak measurement [19]. Many details of a 
practical parity violating experiment to measure neutron densities have been discussed in a long 
paper [10]. 

The neutron radius of 208 Pb, R n , has important implications for astrophysics. There is a strong 
correlation between R n and the pressure of neutron matter P at densities near 0. 1 fm~ 3 (about 2/3 
of nuclear density) [21]. A larger P will push neutrons out against surface tension and increase R n . 
Therefore measuring R n constrains the equation of state (EOS) — pressure as a function of density 
— of neutron matter. 

Recently Hebeler et al. [22] used chiral perturbation theory to calculate the EOS of neutron 
matter including important contributions from very interesting three neutron forces. We have some 
information on isospin 1/2 three nucleon forces from mass 3 nuclei ( 3 He, 3 H) and proton-deuteron 
scattering. However, our experimental information on three neutron forces is limited. From their 
EOS, they predict R n - R p = 0. 17 ± 0.03 fm. Here R p is the known proton radius of 208 Pb. Monte 
Carlo calculations by Carlson et al. also find sensitivity to three neutron forces [23]. Therefore, 
measuring R n provides an important check of fundamental neutron matter calculations, and con- 
strains three neutron forces. 

The correlation between R n and the radius of a neutron star r^s is also very interesting [24]. 
In general, a larger R n implies a stiff er EOS, with a larger pressure, that will also suggest r^s is 
larger. Note that this correlation is between objects that differ in size by 18 orders of magnitude 
from R n pa 5.5 fm to r^s ~ 10 km. We discuss observations of r#s in Section 3. 

The EOS of neutron matter is closely related to the symmetry energy S. This describes how the 
energy of nuclear matter rises as one goes away from equal numbers of neutrons and protons. There 
is a strong correlation between R n and the density dependence of the symmetry energy dS/dn, 
with n the baryon density. The symmetry energy can be probed in heavy ion collisions [25] . For 
example, dS/dn has been extracted from isospin diffusion data [26] using a transport model. 

The symmetry energy 5 helps determine the composition of a neutron star. A large S, at high 
density, implies a large proton fraction Y p that will allow the direct URCA process of rapid neutrino 
cooling. If R n —R p is large, it is likely that massive neutron stars will cool quickly by direct URCA 
[27]. In addition, the transition density from solid neutron star crust to the liquid interior is strongly 
correlated with R n —R p [28]. 

Finally, atomic parity violation (APV) is sensitive to R n [29], [30], [10]. Parity violation in- 
volves the overlap of atomic electrons with the weak charge of the nucleus, and this is primarily 
earned by the neutrons. Furthermore, because of relativistic effects the electronic wave function 
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can vary rapidly over the nucleus. Therefore, the APV signal depends on where the neutrons are 
and on R n . A future low energy test of the standard model may involve the combination of a precise 
APV experiment along with PV electron scattering to constrain R n . Alternatively, measuring APV 
for a range of isotopes can provide information on neutron densities [31]. 

2.2 The Lead Radius Experiment (PREX) 

We now discuss a direct measurement of R n . Parity violation provides a model independent 
probe of neutrons, because the Z° boson couples to the weak charge, and the weak charge of a 
proton <2{y = 1 —4 sin 2 ©iv ~ 0.05 is much smaller than the weak charge of a neutron Qfa = — 1. 
Here &w is the weak mixing angle. 

The Lead Radius Experiment (PREX) at Jefferson Laboratory [9] measures the parity violating 
asymmetry A pv for elastic electron scattering from 208 Pb. The asymmetry A pv is the fractional cross 
section difference for scattering positive (+), or negative (-), helicity electrons, 

da I _ da I 

a dQ. I + dQ. I ~ / ry i \ 

A F* - dai 7¥| • K ' 

da\+^~ dQi- 

In Born approximation, A pv arrises from the interference of a weak amplitude of order the Fermi 
constant Gp , and an electromagnetic amplitude of order the fine structure constant a over the 
square of the momentum transfer q 2 [19], 
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Here the weak form factor Fw(q 2 ) is the Fourier transform of the weak charge density pw{r), that 
is essentially the neutron density, Fw(q 2 ) = J d 3 r sm ^ pw(r), see Fig. 2. 

Likewise, the electromagnetic form factor F c h{q 2 ) is the Fourier transform of the (electro- 
magnetic) charge density p c h( r )- This is known from elastic electron scattering [11]. Therefore, 
measuring A pv as a function of q allows one to map out the neutron density p„(r). Note that, for 
a heavy nucleus, there are important corrections to Eq. 2.2 from Coulomb distortions. However, 
these have been calculated exactly by solving the Dirac equation for an electron moving in both a 
Coulomb potential of order 25 MeV and a weak axial vector potential of order electron volts [32]. 
Therefore, even with Coulomb distortions, one can accurately determine neutron densities. Note 
that this purely electroweak reaction is free from most strong interaction uncertainties. 

The PREX experiment measures A pv for 1.05 GeV electrons elastically scattered from 208 Pb at 
laboratory angles near five degrees. The first measurement yielded A pv = 0.656 ±0.060 (statisti- 
cal) ±0.014 (systematic) ppm [9]. From this the rms neutron radius R n minus proton radius R p for 
208 Pb was found to be R n —R p = 0.33lQ'}g fm. See also ref. [33] for more details of this analyis. 
A second PREX run is now approved to accumulate more statistics and reach the original goal of 
determining R n to 1% (±0.05 fm). 

In addition to PREX, many other parity violating measurements of neutron densities are pos- 
sible, see for example [34]. Measuring R n in 48 Ca is particularly attractive. First, 48 Ca has a higher 
experimental figure of merit than 208 Pb. Therefore a 48 Ca measurement may take less beam time 
than for 208 Pb. Not only does 48 Ca have a large neutron excess, it is also relatively light. With only 
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Figure 2: Densities of 208 Pb in a relativistic mean field model [20]. The lower dashed (red) curve shows 
the point proton density while the upper dashed (black) curve is the neutron density. The (electromagnetic) 
charge density is the lower solid (red) curve while the upper solid (black) curve shows the weak charge 
density that is measured in the PREX experiment. 

48 nucleons, microscopic coupled cluster calculations [35], or no core shell model calculations 
[36], may be feasible for 48 Ca that are presently not feasible for 2()8 Pb. Note that these microscopic 
calculations may have important contributions from three nucleon forces. This will allow one to 
make microscopic predictions for the neutron density and relate a measured R n to three nucleon 
forces and in particular to very interesting three neutron forces. 

3. Electromagnetic observations of neutron star radii 

Neutron stars are collapsed stellar objects that are formed in supernova explosions. They are 
more massive than the sun but have radii of order 10 km. This makes them an extraordinary 18 
orders of magnitude larger than a 208 Pb nucleus and 55 orders of magnitude more massive. There 
is a very interesting relationship between the neutron radius of 208 Pb, of order 6 fm, and the radius 
of a neutron star, of order 10 km [37]. This involves a breathtaking extrapolation of 18 orders 
of magnitude in size, or 55 orders of magnitude in mass. Nevertheless, both in the laboratory 
and in astrophysics, it is the same neutrons, the same strong interactions, the same neutron rich 
matter, and the same equation of state. A measurement in one domain, be it astrophysics or in 
the laboratory, can have important implications in the other domain. This is the real strength of a 
multi-messenger approach to neutron rich matter. 

3.1 Neutron star radii 

The structure of a neutron star can be calculated with the Tolman-Oppenheimer-Volkoff Equa- 
tions of General Relativity [38] and is completely determined by the equation of state of neutron 
rich matter. The radius of a neutron star depends on the pressure of neutron matter at normal nu- 
clear density and above, because the central density of a neutron star can be a few or more times 
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that of normal nuclear density. A higher pressure will lead to a larger radius. It is important to have 
both low density information on the equation of state from PREX, and high density information 
from measurements of neutron star radii. This can constrain any possible density dependence of 
the equation of state from an interesting phase transition to a possible high density exotic phase 
such as quark matter, strange matter, or a color superconductor. For example, if the 208 Pb radius is 
relatively large, this shows the EOS is stiff at low density (has a high pressure). If at the same time, 
neutron stars have relatively small radii, than the high density EOS is soft with a low pressure. This 
softening of the EOS with density could strongly suggest a phase transition to a soft high density 
exotic phase. 

The radius of a neutron star r^s can be deduced from X-ray measurements of luminosity L and 
surface temperature T , 

L = AKr 2 NS a SB T\ (3.1) 

with Osb the Stefan Boltzmann constant. This allows one to deduce the surface area of a neutron 
star. However, there are important complications. First, one needs an accurate distance to the neu- 
tron star. Second, Eq. 3. 1 assumes a black-body and there are important non-blackbody collections 
that must be determined from models of neutron star atmospheres. Finally, gravity is so strong that 
space near a neutron star is strongly curved. If one looks at the front of a neutron star, one will also 
see about 30% of the back, because of the curvature of space. Therefore the surface area that one 
observes in Eq. 3. 1 depends on the curvature of space and the mass of the star. Thus what started 
out as a measurement of just the radius is, in fact, a measurement of a combination of mass and 
radius. 

Recently Steiner, Lattimer, and Brown have deduced masses and radii [39] from combined 
observations of six neutron stars in two classes: 1) X-ray bursts, and 2) neutron stars in globular 
clusters. They conclude that observations favor a stiff high density equation of state that can support 
neutron stars with a maximum mass near 2 M Q and that the equation of state is soft at low densities 
so that a 1.4 M Q neutron star has a radius near 12 km. They go on to predict that the neutron 
minus proton root mean square radius in 208 Pb will be R n — R p = 0.15 ± 0.02 fm. Note that this 
is a prediction for a nucleus based on an equation of state deduced from X-ray observations of 
neutron stars. The Steiner et al. paper [39] is potentially controversial because their results depend 
on, among other things, the model assumed for X-ray bursts. Ozel et al. use a different model for 
X-ray bursts, and get very small neutron star radii near 10 km or below [40]. The high density EOS 
that Ozel et al. deduce from these observations is significantly softer than Steiner et al's, suggesting 
a phase transition to an exotic phase [41]. Thus, the difference between Ozel et al's result of 10 
km, that likely implies a phase transition, and Steiner et al's result of 12 km, that does not suggest a 
phase transition, is very significant. Clearly these observations of neutron star radii have potentially 
very important implications for the properties of neutron rich matter and the high density phases of 
QCD. 

3.2 Solid neutron rich matter and neutron star crust cooling 

We go on to discuss the neutron star crust and crust cooling. The neutron radius of 208 Pb has 
implications for neutron star structure, in addition to the star's radius. Neutron stars have solid 
crusts over liquid cores, while a heavy nucleus is expected to have a neutron rich skin. Neutron 
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star crust consists of a relativistic fermi gas of electrons, a crystal lattice of neutron rich ions, and, 
in general, a neutron gas. For a review see [42]. Both the solid crust of the star, and the skin of 
the nucleus, are made of neutron rich matter at similar, slightly subnuclear, densities. The common 
unknown is the equation of state of neutron matter. A thick neutron skin in 208 Pb, means a high 
pressure where the energy rises rapidly with density. This quickly favors the transition to a uniform 
liquid phase. Therefore, we find a strong correlation between the neutron skin thickness, measured 
by PREX, and the transition density in neutron stars from solid crust to liquid interior [28]. 

Perhaps the presence of a solid crust deserves comment. We expect solids to form in compact 
stars, both white dwarfs and neutron stars. This may sound surprising since stars are often com- 
posed of not liquids or gasses, but plasmas. Nevertheless, the plasmas can be so dense that the ions 
actually freeze. Recently, there are significant new observations of how white dwarfs freeze [43] 
and of how the solid crust of accreting neutron stars cools [44, 45, 46]. In addition, crystalization 
has been observed in complex laboratory plasmas, see for example [47, 48]. We hope to extract 
important new information, on the structure of solids in stars, by comparing these observations to 
large scale molecular dynamics simulations of freezing in both white dwarfs [49] and neutron stars 
[50]. 

4. Gravitational Waves 

We turn now to gravitational wave observations of neutron rich matter. Albert Einstein, almost 
100 years ago, predicted the oscillation of space and time known as gravitational waves (GW). 
Within a few years, with the operation of Advanced LIGO [51], Advanced VIRGO [52] and other 
sensitive interferometers, we anticipate the historic detection of GW. This will be a remarkable 
achievement and open a new window on the universe and on neutron rich matter. 

The first GW that are detected will likely come from the merger of two neutron stars. The 
rate of such mergers can be estimated from known binary systems [53]. During a merger the GW 
signal has a so called chirp form where the frequency rises as the two neutron stars spiral closer 
together. Deviations of this wave form from that expected for two point masses may allow one to 
deduce the equation of state of neutron rich matter and measure the radius of a neutron star rws 
[54]. Alternatively one may be able to observe the frequency of oscillations of the hyper-massive 
neutron star just before it collapses to a black hole. This frequency depends on the radius of the 
maximum mass neutron star [55]. However, either approach may require high signal to noise data 
from relatively nearby mergers. Note that this information on r^s from GW is independent of the 
X-ray observations discussed in Sec. 3 and of any systematic errors associated with the modeling 
of X-ray bursts. 

Continuous GW signals can also be detected, see for example [56]. Indeed Bildstein and others 
[57] have speculated that some neutron stars in binary systems may radiate angular momentum in 
continuous GW at the same rate that angular momentum is gained from accretion. This would 
explain why the fastest observed neutron stars are only spinning at about half of the breakup rate. 
There are several very active ongoing and near future searches for continuous gravitational waves 
at LIGO, VIRGO and other detectors, see for example [58]. No signal has yet been detected. 
However, sensitive upper limits have been set. These limits constrain the shape of neutron stars. 
In some cases the star's elipticity £, which is that fractional difference in moments of inertia e = 
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{h —h)/h '^ observed to be less than a part per million or even smaller. Here h,h, and It, are the 
principle moments of inertia. 

In general, the amplitude of any continuous signal is much weaker than a burst signal. How- 
ever, one can gain sensitivity to a continuous signal by coherently (or semi-coherently) integrating 
over a large observation time, see for example [59]. Note that searches for continuous GW can 
be very computationally intensive because one must search over an extremely large space of pa- 
rameters that may include the source frequency, how that frequency changes with time, the source 
location on the sky, etc. The Einstein at home distributed computing project uses spare cycles on 
the computers of a large number of volunteers to search for continuous GW [60]. 

Strong GW sources often involve large accelerations of large amounts of neutron rich matter. 
Indeed the requirements for a strong source of continuous GW, at LIGO frequencies, places ex- 
traordinary demands on neutron rich matter. Generating GW sounds easy. Place a mass on a stick 
and shake vigorously. However to have a detectable source, one may need not only a large mass, 
but also a very strong stick. The stick is needed to help produce large accelerations. Since others 
have discussed large masses, let us focus here on the strong stick. 

An asymmetric mass on a rapidly rotating neutron star produces a time dependent mass 
quadrupole moment that radiates gravitational waves. However, one needs a way (strong stick) 
to hold the mass up. Magnetic fields can support mountains, see for example [61]. However, it may 
require large internal magnetic fields. Furthermore, if a star also has a large external dipole field, 
electromagnetic radiation may rapidly spin the star down and reduce the GW radiation. 

Alternatively, mountains can be supported by the solid neutron star crust. Recently we per- 
formed large scale MD simulations of the strength of neutron star crust [62, 63]. A strong crust can 
support large deformations or "mountains" on neutron stars, see also [64], that will radiate strong 
GW. How large can a neutron star mountain be before it collapses under the extreme gravity? This 
depends on the strength of the crust. We performed large scale MD simulations of crust breaking, 
where a sample was strained by moving top and bottom layers of frozen ions in opposite directions 
[62]. These simulations involve up to 12 million ions and explore the effects of defects, impurities, 
and grain boundaries on the breaking stress. 

We find that neutron star crust is very strong because the high pressure prevents the forma- 
tion of voids or fractures and because the long range coulomb interactions insure many redundant 
"bounds" between planes of ions. Neutron star crust is the strongest material known, according 
to our simulations. The breaking stress is 10 billion times larger than that for steel. This is very 
promising for GW searches because it shows that large mountains are possible, and these could 
produce detectable signals. 

Continuous GW can also be produced by r-mode oscillations of a rotating neutron star [65]. 
Consider a surface wave on a rapidly rotating neutron star that is moving slowly in a direction 
opposite to the stars rotation. This wave, in the laboratory frame, will appear to be moving in the 
direction of the star's rotation. The back reaction force on the wave from GW radiation will always 
act to slow the wave in the laboratory frame. However, in this case slowing in the lab frame will 
actually speed up the wave in the rotating star's frame and increase its amplitude. Thus the wave 
can be unstable with respect to gravitational wave radiation. The r-modes are collective oscillations 
on rotating neutron stars that can also be unstable to GW radiation [65]. If an r-mode is unstable the 
amplitude will grow large and rotational kinetic energy can be radiated away as GW. This will slow 
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the rotation rate. Note, that the physics of large amplitude r-mode oscillations can be complicated, 
see for example [66, 67]. 

The stability of r-modes depends on the amount of dissipation from, for example, the bulk and 
shear viscosities of neutron rich matter. If dissipation is large then the amplitude of the r-modes 
will stay small and rapid rotation of a neutron star is possible. Alternatively, if dissipation is small 
then the r-modes may be unstable and GW radiation from the modes may limit the rotation rate 
of the star. Unfortunately the stability of r-modes has proved to be a complex subject that may be 
sensitive to subtle dissipation properties of neutron rich matter, be it in a nucleon phase [68, 69] or 
in more exotic quark and gluon phases [70, 71, 72]. 

We give one example of a possible source of dissipation for the r-modes. The shear viscosity of 
conventional complex fluids, with large non-spherical molecules, can be orders of magnitude larger 
than that for normal fluids. This suggests that the shear viscosity of nuclear pasta, with long rod 
like shapes as seen in Fig. 1, could be large. In ref. [73], we have calculated the shear viscosity of 
nuclear pasta using large scale molecular dynamics simulations. The shear viscosity is dominated 
by momentum carried by electrons and although the electron mean free path is determined by 
electron-pasta scattering, we find no dramatic differences from a conventional phase with spherical 
nuclei. Therefore, we find that the shear viscosity of nuclear pasta is not very different from that 
for more conventional matter with nearly spherical nuclei. However there could be other sources 
of dissipation. We do not yet have a complete understanding of when the r-modes may be stable 
and when they are unstable. 

To conclude this section, there is a great deal of interest in gravitational waves (GW) from 
neutron stars and there are many ongoing searches. One is interested in both burst sources, for 
example from neutron star mergers, and continuous sources from mountains or collective modes. 
Gravitational wave radiation depends on the equation of state of neutron rich matter. In addition, it 
can also depend on other more detailed properties including the breaking strain of solid phases and 
the bulk and shear viscosities. 

5. Supernova neutrinos and neutron rich matter 

Neutrinos provide yet another, qualitatively different, probe of neutron rich matter. Core col- 
lapse supernovae (SN) are gigantic stellar explosions that convert as much as 0.2M© of mass into 
10 58 neutrinos [74]. We detected about 20 neutrinos from SN1987a [75, 76]. The next galactic SN 
should be very exciting. There are many new underground experiments to search for dark matter, 
double beta decay, solar neutrinos, proton decay, oscillation of accelerator neutrinos, etc. Some of 
these experiments can be very sensitive to SN neutrinos. For example, the coming ton scale dark 
matter experiments can be sensitive to SN neutrinos via neutrino-nucleus elastic scattering [77]. 

In a SN, neutrinos are emitted from the surface of last scattering known as the neutrino-sphere. 
At the neutrino-sphere, the neutrino mean free path is comparable to the size of the system. The 
temperature is of order 4 MeV, based on the energies of the SN1987a events. The density is w 10 11 
g/cm 3 , 1/1000 to 1/100 of normal nuclear density. This follows from the mean free path implied 
by known neutrino cross sections. Thus, the neutrino-sphere is a low density warm gas of neutron 
rich matter. The pressure, composition, and long wave-length neutrino response [78] of this region 
can be calculated from a model independent Virial expansion [79, 80]. This expansion is based on 
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nucleon-nucleon, nucleon-alpha, and alpha-alpha elastic scattering phase shifts. Neutrino-sphere 
conditions can be simulated in the laboratory with heavy ion (HI) collisions. For example during a 
peripheral collision, fragments emitted with velocities intermediate between those of the projectile 
and target are coming from a warm low density region, see for example [81]. 

The next galactic SN should provide a great deal of information on neutrino properties such 
as oscillations, masses, and mixing angles. We give one example of an important SN neutrino 
observable related to nucleosynthesis. About half of the elements heavier than 56 Fe, including 
gold and uranium, are thought to be made in the r-process of rapid neutron capture nucleosynthesis 
[83]. Here seed nuclei rapidly capture many neutrons to produce very neutron rich nuclei that then 
beta decay several times to produce heavy elements such as gold. At this time, the preferred site 
for the r-process is the neutrino driven wind during a supernova, see for example [84] . Here some 
baryons are blown off of the proto-neutron star by the intense neutrino flux. Nucleosynthesis in this 
wind depends on the entropy, the expansion time scale, and most importantly the proton fraction 
Y p (number of protons divided by the total number of protons and neutrons). The proton fraction is 
set by the relative rates of neutrino and antineutrino capture. 

v e + n— >p + e (5.1) 

v e + p->n + e + (5.2) 

The cross sections for these reactions depend on neutrino and antineutrino energies. Therefore 
one should measure the difference AE between the average energy of electron antineutrinos and 
neutrinos, 

AE = (Ey e ) - (E Ve ) . (5.3) 

If AE is large, the wind can be neutron rich. However if AE is small, the wind will be proton 
rich and this is likely a serious problem for r-process nucleosynthesis in the wind [85]. The wind 
has been the preferred r-process site because, for example, it is a known explosive environment 
that occurs relatively often so it can easily supply enough r-process material. However, present 
simulations of the wind do not have a large enough ratio of free neutrons to seed nuclei. One can 
increase this ratio by lowering Y p (and increasing the number of neutrons) or by increasing the 
entropy (and destroying some of the seed nuclei). 

A proton in neutron rich matter is more tightly bound than a neutron. This is because of the 
symmetry energy and this binding energy difference can increase the cross section for Eq. 5.1 and 
decrease the cross section for Eq. 5.2, [82]. This could change the Y p in the neutrino driven wind 
and impact nucleosynthesis. 

Finally, calculations of r-process nucleosynthesis yields also require important nuclear struc- 
ture input including masses of neutron rich nuclei, beta decay half lives, and neutron capture cross 
sections [83]. New radioactive beam accelerators, such as FRIB, can directly produce some of the 
neutron rich nuclei involved in the r-process and provide important nuclear structure data. However 
this nuclear data may not, by itself, directly address important astrophysical questions involving the 
site of the r-process and the source of the neutrons. 

In any case, measuring AE, during the next galactic SN, should provide an important diagnos- 
tic on conditions in the neutrino driven wind. At present the site of the r-process is unknown [86]. 
One alternative site for the r-process is tidally ejected material during neutron star mergers, see for 
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example [87] [88] [89]. If this is the site, there may be an observable electromagnetic signal associ- 
ated with the radioactive decay of r-process nuclei, see for example [90]. The Spaniards started a 
quest to find the city of Eldorado. Perhaps with SN neutrinos or electromagnetic observations of 
neutron star mergers, we can end this quest and find the source of the chemical element gold. 

6. Conclusions: neutron rich matter 

Neutron rich matter is at the heart of many fundamental questions in Nuclear Physics and 
Astrophysics. What are the high density phases of QCD? Where did the chemical elements come 
from? What is the structure of many compact and energetic objects in the heavens, and what 
determines their electromagnetic, neutrino, and gravitational- wave radiations? Moreover, neutron 
rich matter is being studied with an extraordinary variety of new tools such as the Facility for Rare 
Isotope Beams (FRIB) and the Laser Interferometer Gravitational Wave Observatory (LIGO). 

We described the Lead Radius Experiment (PREX) that uses parity violating electron scatter- 
ing to measure the neutron radius in 208 Pb. This has important implications for neutron stars and 
their crusts. We discussed X-ray observations of neutron star radii that also have important implica- 
tions for neutron rich matter. Gravitational waves (GW) from sources such as neutron star mergers 
and rotating neutron star mountains open a new window on neutron rich matter. Using large scale 
molecular dynamics simulations, we found neutron star crust to be the strongest material known, 
some 10 billion times stronger than steel. It can support mountains on rotating neutron stars large 
enough to generate detectable gravitational waves. Finally, neutrinos from core collapse super- 
novae (SN) provide another, qualitatively different probe of neutron rich matter. Neutrinos come 
from the neutrino-sphere that is a low density warm gas phase of neutron rich matter. Neutrino- 
sphere conditions can be simulated in the laboratory with heavy ion collisions. Observations of 
neutrinos can probe nucleosyntheses in SN while SN simulations depend on the equation of state 
(EOS) of neutron rich matter. 

In conclusion, multi-messenger astronomy is based on the widely held belief that combining 
astronomical observations using photons, gravitational waves, and neutrinos will fundamentally 
advance our knowledge of compact and energetic objects in the heavens. Compact objects such 
as neutron stars are, in fact, giant nuclei, even if they are an extraordinary 18 orders of magnitude 
larger than a 208 Pb nucleus. Nevertheless, both in the laboratory and in Astrophysics, these ob- 
jects are made of the same neutrons, that undergo the same strong interactions, and have the same 
equation of state. A measurement in one domain, be it Astrophysics or the laboratory, can have im- 
portant implications in the other domain. Therefore we can generalize multi-messenger astronomy 
to multi-messenger observations of neutron rich matter. We believe that combing astronomical 
observations using photons, GW, and neutrinos, with laboratory experiments on nuclei, heavy ion 
collisions, and radioactive beams will fundamentally advance our knowledge of the heavens, the 
dense phases of QCD, the origin of the elements, and of neutron rich matter. 
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